Antiferromagnetic and structural transitions in the superoxide KO2 from first 
principles: A 2p-electron system with spin-orbital-lattice coupling 
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KO2 exhibits concomitant antiferromagnetic (AFM) and structural transitions, both of which 
originate from the open-shell 2p electrons of molecules. The structural transition is accompanied 
by the coherent tilting of O^" molecular axes. The interplay among the spin-orbital-lattice degrees 
of freedom in KO2 is investigated by employing the first-principles electronic structure theory and 
the kinetic-exchange interaction scheme. We have shown that the insulating nature of the high 
symmetry phase of KO2 at high temperature (T) arises from the combined effect of the spin-orbit 
coupling and the strong Coulomb correlation of O 2p electrons. In contrast, for the low symmetry 
phase of KO2 at low T with the tilted molecular axes, the band gap and the orbital ordering 
are driven by the combined effects of the crystal-field and the strong Coulomb correlation. We have 
verified that the emergence of the O 2p ferro-orbital ordering is essential to achieve the observed 
AFM structure for K0 2 . 

PACS numbers: 75.50.Ee, 71.70.Ej, 71.15.Mb 
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Magnetism due to the correlated 2p electrons has at- 
tracted revived attention for the possibility of the new 
kinds of the magnetic informative materials.—^ Some 2p 
magnetic oxides exhibit the structural phase transition 
concomitantly with the magnetic phase transitional^ 
Solid oxygen is a typical example, which has both an- 
tiferromagnetic (AFM) and structural phase transitions 
below 24 K4 Alkali superoxides, A0 2 (^4=Na, K, Rb), 
which are of our present interest, belong to another 
example.— £ In AO2, one alkali-metal atom provides one 
electron to an oxygen molecule, and thereby each 
anion has nine electrons at the 2p molecular levels with 
the electronic configuration of a^it^n^ (see the inset of 
Fig. [2]) The partially occupied antibonding ir g molecu- 
lar states play the most important role in determining the 
electronic and magnetic properties of alkali superoxides. 8 
One hole in n fl generates the magnetic moment of 1 fi b 
for each O^ . The degeneracy of the Tr g level is expected 
to be lifted by lowering the crystal symmetry, as occurs 
due to the Jahn- Teller effect 4 In fact, it was suggested 
that, for KO2, the symmetry lowering would occur via 
coherent tilting of the molecular axes, the so-called 
magnetogyration, which invokes the accompanying AFM 
ordering 4££ 

At room temperature, KO2 crystallizes in the tetrag- 
onal structure of CaC2 type, in which the molec- 
ular bond axes are parallel to the z-axis (Fig. [T^a)). 4 
KO2 retains this structure down to 197 K and exhibits 
the paramagnetic behavior. Upon cooling, O^ molecu- 
lar bond axes seem to tilt uniformly by ~ 20° to have 
a lower crystal (monoclinic) symmetry. The magnetic 
phase is still paramagnetic down to 7 K. Below 7 K, the 
AFM ordering emerges in the triclinic crystal structure 
with the uniform tilting of 0^~ molecular bond axes by 
~ 30° A According to neutron experiment^ the AFM 
phase has the magnetic structure having opposite spin 
arrangements along the ^-direction between two oxygen 
layers of 01 and 02. This feature in KO2 reflects the 



strong interplay among spin, orbital, and lattice degrees 
of freedom, as in rare-earth manganites^ 

There have been several theoretical reports to study 
the coupled structural and magnetic transitions in 
K0 2 ,-> 8 -&±^ but those studies were mostly qualitative 
and lacked the quantitative description of the electronic 
structures for the low symmetry phase of KO2. Even 
the direction of the tilted molecular bond axis is still un- 
certain between two possibilities. The first one is [Rl] 
in Fig. [ljb), in which the molecular bonds are rotated 
around the [100] axis. The second one is [R2] in Fig.QJc), 
in which the molecular bonds are rotated around [110] 
axis. Moreover, the recent band structure study for the 
high symmetry phase of KO2 ((Fig. Qla)) in the local 
density approximation (LDA)^ reveals that the degener- 
ate TT g states do not split, resulting in the metallic nature, 
which is contradictory to the insulating nature of KO2. 

In this Rapid Communication, we have investigated 
comprehensively the electronic structures of both high 
and low symmetry phases of KO2, and explored the 




FIG. 1: (Color online) (a) The tetragonal [T] structure of KO2 
with unrotated molecular bond axes. The molecular axes 
are along the z-direction. K atoms are in blue and O atoms in 
red. (b) [Rl] with molecular bonds rotated by ~ 30° around 
the [100] (x') axis with bond axes along the z'-direction. (c) 
[R2] with molecular bonds rotated by ~ 30° around the [110] 
(as") axis with bond axes along the z"-direction. There are 
two independent types of oxygen 01 and 02 for each struc- 
ture. 
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FIG. 2: (Color online) The DOS in the GGA for the FM phase 
of KO2 with [T] structure. Dotted (blue) and solid (red) lines 
represent the total DOS and the partial DOS of oxygen 2p 
states, respectively. The upper and the lower panels corre- 
spond to the majority and minority spin DOSs. Inset shows 
the 2p molecular levels of OJ. 



FIG. 3: (Color online) The PDOS of oxygen 2p states for the 
FM phase of KO2 with [T] structure. The dotted (blue) and 
solid (red) lines represent the DOSs in the GGA+C/ and the 
GGA+C/+SO, respectively ((7=6.53 eV). The tt 9 states split 
into 7r OT= i,7r m= _i in the GGA+(7+SO. Charge densities of 
7r m =i states are plotted on the xy (001) (the upper panel) 
and on the xz (010) plane (the lower panel). 



mechanism of the interplay between the spin and lattice 
degrees of freedom in correlated 2p electron systems. Fur- 
ther, we have examined the origins of the band gap open- 
ing and the exchange interactions between 0.7 molecules 
in K0 2 . 

We have employed the full-potential augmented plane 
wave (FLAPW) band method 13 implemented in WIEN2k 
packaged For the exchange-correlation potential, the 
generalized gradient approximation (GGA)i£ was used. 
We also incorporate the on-site Coulomb interaction U 
between the oxygen 2p electrons (GGA+C7) and the 
spin-orbit (SO) effect as a second variational procedure 
(GGA+C7+SO). The valence wave functions inside the 
muffin-tin spheres are expanded in spherical harmon- 
ics up to Z m ax=10, and the wave function in the in- 
terstitial region is expanded with plane waves up to 
K ma x = 7/i?MT, where Rmt is the smallest muffin-tin 
sphere radius. -Rmt's were set as 2.00 (a.u.) for K and 
1.23 (a.u.) for Oxygen, respectively. The charge density 
was expanded with plane waves up to G max =12 (a.u.) -1 . 
We have used 200 k-points inside the first Brillouin zone. 
For [Rl] and [R2] in Fig.[TJ the 07 molecular bonds were 
assumed to be rotated uniformly by 30° from [T] struc- 
ture around the [100] and [110] axis, respectively. The 
experimental lattice parameters a=4.030A and c=6.697A 
and the bond length d oo =1.306A were usedJ£ The spin 
direction was chosen to be perpendicular to the molecular 
axis following the experiment^ 

Figure [5] provides the density of states (DOS) in the 
GGA band calculation for the FM phase of KO2 with the 
high symmetry un-rotated structure [T] . Oxygen 2p par- 
tial DOS (PDOS) shows clear molecular level splittings 
among <7 9 ,7r u ,7r s and a u in agreement with literature^ 
The finite DOS at the Fermi level Ef in the minority spin 
■Kg states produces the half- metallic character for KO2, 



which reflects the failure of the GGA. The degenerate tt 9 
states, which are composed of anti-bonding n x and tt v 
states, should split for an insulating phase of K02- To 
understand the splitting mechanism of the 7r ff states, we 
have considered the effects of the SO and the Coulomb 
correlation of O 2p electrons. The importance of the 
Coulomb correlation effect of O 2p electrons was exam- 
ined recently in a similar 2p magnetic oxide Rb40 6 i 17 i 18 

Figure 1 shows the DOS both in the GGA+C7 and 
GGA+C7+SO for the FM phase of K0 2 with [T] struc- 
ture. Note that the ir g states near Ep- do not split in the 
GGA+C7 with [7=6.53 eV. They become split only when 
the SO is included in the GGA+C7: ir g split into ir m= i 
and Tr m= _i which are mixed states of ir x and n y . In the 
right of Fig. [21 the charge densities for the occupied 7r m= i 
states are plotted on the (001) and (010) planes. These 
charge densities manifest the azimuthally symmetric na- 
ture of 7r m= i states, which demonstrates that the split- 
ting of Wg states originates from the SO effect. In fact, 
the large SO effect results from the Coulomb correlation, 
which localizes the 2p electrons to generate the substan- 
tial orbital magnetic moment in the ir g states. The mag- 
nitude of the orbital magnetic moment is as much as 
0.570 /is per 07 molecule, which is comparable to that 
of the spin magnetic moment of 1.0 /is- The direction 
of the orbital magnetic moment turns out to be almost 
parallel to the molecular axis, i.e., perpendicular to that 
of the spin magnetic moment. We have checked that the 
electronic structure of AFM KO2 with [T] structure is 
close to that of FM KO2, and so the band gap at E^ 
opens in the same way. 

For KO2 with [Rl] and [R2] structures, the crystal field 
effect will be activated due to tilting of O^ molecular 
axes toward K + . Figure IDJa) shows the local PDOSs 
and the spin densities in the GGA+C7+SO for the AFM 
phase of KO2 with [Rl] structure. Notable feature in 
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FIG. 5: (Color online) (a) The in-plane exchange constants, 
Ji and J2, on the rsy-plane for KO2 with [R2], where the un- 
occupied tt x >i and occupied ir y » states are represented by red 
and green, respectively, (b) The inter-plane exchange con- 
stant J3 on the (110) plane. The paths A and B correspond 
to the direct and indirect (via K + ) hopping channels, respec- 
tively. 



FIG. 4: (Color online) (a) The local PDOS of the O 2p states 
in the GGA+C/+SO ((7=6.53 eV) for the AFM phase of K0 2 
with [Rl] structure. The spin densities of unoccupied n y i 
states are plotted on the (100) and (200) planes. Blue and red 
colors in the spin densities represent opposite spins, (b) The 
local PDOS of the O 2p states in the GGA+I7+SO ((7=6.53 
eV) for the AFM phase of K0 2 with [R2] structure. The spin 
density of occupied minority spin 7iy/ states is plotted on the 
(110) plane. Black solid lines in the DOS plots represent the 
PDOS of K 4p states. 



the DOS is the opening of band gap at Ej?. The GGA 
band calculation does not produce the band gap, whereas 
both the GGA+C7 and the GGA+C/+SO yield almost 
the same insulating band structure. This implies that, 
in producing the band gap, (i) the symmetry lowering 
just by tilting of molecular axes is not sufficient, (ii) the 
Coulomb correlation effect is essential, and (iii) the SO 
effect for [Rl] is not as important as that for [T] . The or- 
bital magnetic moment is almost quenched having only 
0.002 [Ib per 0^~ molecule. Therefore, it is the combined 
effect of the Coulomb U and the crystal field from K + 
cations that splits the degenerate 7r s states into the di- 
rectional orbitals, ir x i and 7r H < states. Interestingly, due 
to the complex interplay of the Coulomb U and the crys- 
tal field effect in the molecular states, the splitting of tt 9 



TABLE I: Magnetic exchange constants [meV] between O^ 
molecules for direct and indirect (via K + ) hopping channels 
in the [Rl] and [R2] structures of KO2. Ji and J2 are the in- 
plane exchange constants along the x and the y-direction, and 
Jz is the inter-plane exchange constant. Positive J represents 
the AFM interaction. 





Ji 


J2 


J3 


Rl (direct) 


0.353 


1.486 


0.209 


Rl (via K+) 


0.023 


0.023 


0.027 


R2 (direct) 


0.419 


0.419 


0.680 


R2 (via K+) 


-0.215 


-0.215 


0.005 



states occurs extraordinarily, i.e., the 7iy-7iv order for 
the majority spin and the opposite 7T x /-7ry order for the 
minority spin states. The half-filled molecular states 7iy 
seem to be strongly affected by Coulomb U, while the 
fully occupied molecular states ir x / seem to be weakly 
affected. The spin densities of unoccupied ir g states plot- 
ted on the (100) and (200) planes demonstrate that they 
really correspond to the ir y > states. The 7ty states in 
each plane exhibit the FM spin and ferro-orbital (FO) 
orderings, while they exhibit the inter-plane AFM spin 
ordering. Namely, the AFM and FO orderings take place 
concurrently. 

Similarly, Fig. BJb) shows the local PDOSs and the 
spin densities in the GGA+(7+SO for AFM K0 2 with 
[R2] structure. For [R2] too, the GGA does not pro- 
duce the band gap, whereas both the GGA+U and the 
GGA+(7+SO produce it. The SO effect is again negli- 
gible so that the orbital magnetic moment is quenched 
having only 0.008 [Ib per 0^~ molecule. Degenerate 7r g 
states are split into tt x " and TT y " due to the combined ef- 
fect of the Coulomb U and the crystal field from K + . The 
spin density of occupied minority spin ir g states plotted 
on the (110) plane demonstrates that they really corre- 
spond to the 7Tj," states. 

Now let us examine the magnetic interaction in KO2 
based on the above electronic structures. The implica- 
tion of Fig. 0] is that, with the tilting of 0^~ molecular 
axes in the AFM [Rl] and [R2] structures, the band gap 
opening and the FO ordering in ir g states occur simul- 
taneously. According to the spin-orbital model for 3d 
transition-metal (TM) oxides^ the FO ordering of TM 
3d states would induce the AFM superexchange inter- 
action between two TM spins, and vice versa. In this 
respect, the results in Fig. 2] look reasonable. The su- 
perexchange, more generally, the kinetic exchange in KO2 
will take place through two channels: the direct hopping 
between two 0^~ molecules and the indirect hopping via 
K + (see Fig. [3]). The unoccupied K 4p states are located 
4 — 6 eV above (see Fig. [5]), and so there is a possibil- 
ity of hopping mechanism through 0^-K + -0^. Indeed, 
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as revealed in Fig. [4] there exists the nonnegligible hy- 
bridization between the O 2p and K Ap states. 

Table U presents the estimated exchange constants 
based on the microscopic calculation of the kinetic ex- 
change interaction for KO2 with the [Rl] and [R2] 
structures^ The 2p electrons in each are assumed 
to have the atomic orbital states as in our band results 
for [Rl] and [R2] . The kinetic exchange interactions both 
from the direct and indirect hoppings were considered in- 
dependently to calculate the exchange constant for each 
channel. For [Rl], the exchange constants are obtained to 
be all AFM, and the dominant channel is the direct hop- 
ping between O^T anions along the y-direction, which is 
AFM (J 2 =1.486 meV). The large J 2 results from the tilt- 
ing of molecular axis along the [010]. The inter-plane 
interactions from both direct and indirect channels yield 
smaller exchange constants (J3=0.209 and 0.027 meV). 
Then the resulting magnetic structure for [Rl] will not be 
consistent with the experimental AFM structure which 
has alternating AFM spins along [001]. 

On the other hand, the kinetic exchange interactions 
for [R2] turn out to be consistent with the experimen- 
tal AFM structure. The dominant one is the inter- 
plane AFM interaction coming from the direct hopping 
(J3=0.680 meV). The in-plane exchange constants are 
0.419 and —0.215 meV for direct and indirect channels, 
respectively. Hence the total exchange constants between 
two nearest neighbor (NN) along x and y-directions 
are 0.204 meV. For each in [R2], there are eight inter- 
plane NN O2 and four in- plane NN 0^~. As a conse- 
quence, much larger J3 than J\ and J 2 will generate the 
interlayer AFM ordering along [001], as is consistent with 
the experimental AFM structure. 

The in-plane FM interactions (—0.215 meV) coming 
from the indirect hopping via K + play an important role 
in stabilizing the in-pane FM ordering. As shown in 



Fig. [5ja), out of two neighboring molecules on the 
xy plane in [R2] , the lobes of tt x " orbital of one molecule 
are toward the intermediate K + , while those of the other 
molecule are away from K + and thereby the lobes of ir y " 
orbital are toward the intermediate K + . This arrange- 
ment of the orbitals suggests that the hopping takes place 
between the occupied 7r y » and unoccupied 7r x » orbitals 
of two neighboring molecules. This is reminiscent 
of the resulting FM superexchange interaction between 
the fully-occupied and half-filled orbitals in TM oxides in 
the framework of the Goodenough-Kanamori- Anderson 
(GKA) rule^i In contrast, for the inter-plane kinetic ex- 
change in Fig. 0Jb), the lobes of ir x » orbitals of both 
molecules are away from K + , and so the indirect kinetic 
exchange results in the weak AFM interaction ( J3=0.005 
meV). 

In conclusion, we have investigated electronic and mag- 
netic structures of KO2 superoxide, the strongly corre- 
lated 2p electron system. We have found that, for the 
correct description of the insulating electronic structure 
of KO2, the SO coupling as well as the large Coulomb cor- 
relation is important for the high symmetry phase, while, 
for the low symmetry phase, the crystal field from K + , 
as well as the large Coulomb correlation, is important. 
The concurrent AFM spin and FO orderings with the 
band-gap opening clearly demonstrate the strong cou- 
pling among the spin-orbital-lattice degrees of freedom 
in KO2. In the low symmetry phase of KO2 with [R2] 
structure, the emergent FO ordering yields the kinetic 
exchange interactions that are consistent with the exper- 
imental AFM structure. 
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